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Abstract 

The triatomic hydrogen ion (H^) has spurred tremendous interest in astrophysics in recent 
decades, and Rydberg states of H3 have also maintained an important role for understanding 

experiments. In a previous study [J. Chem. Phys. 133, 234302 (2010)], radiative transi- 
tions between neutral H3 Rydberg states were calculated at wavelengths near 7 microns, and 
could be compared with mid-infrared laser lines observed in hydrogen/rare gas discharges. The 
present study extends the investigation to wavelengths near 10-13 microns. Rydberg states 
of D3 are also treated. 

Introduction 

Although triatomic hydrogen (H3) and its ion (H^) are the simplest polyatomic molecules, they 
have continued to attract intense interest in diverse contexts, ranging from chemistry to astronomy, 
ever since their discovery. plays an important role in astrophysics since it acts as a proton 
donor in chemical reactions occurring in interstellar clouds.J-^ Furthermore, this ion also helps 
to characterize Jupiter's atmosphere from afar.-"^ is the dominant positively charged ion in 
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molecular hydrogen plasmas and was first identified in 191 1 by J. J. Thomson with an early form of 
mass spectrometry.- Without a stable electronic excited state and a permanent dipole moment, 
cannot be observed by electronic spectroscopy or rotational spectroscopy. Therefore, an infrared 
rotation- vibration spectrum is the only means to observe this ion. The first observation was carried 
out by T. Oka in 1980.-^ By 2012, more than 600 low-lying ro vibrational states of had been 
identified. The good agreement achieved between the experimental spectrum and a first-principles 
calculation provided a benchmark for calculations on other polyatomic molecules such as water. 

One of the biggest surprises among the properties of this simple ion is its dissociative 
recombination (DR) rate, which is important for understanding observations of in diffuse in- 
terstellar clouds.- Until 2003, the DR process, H^-i-e"— >H3— )■ H2-1-H or H-i-H-i-H, has been studied 
in several different experiments, and had an order of magnitude discrepancy with theoretical ex- 
pectation at that time. Building on the previous work of Schneider, Orel, and Suzor-Weiner, - 
Kokoouline and Greene showed'^'^ that intermediate Rydberg states of H3 play an important role 
in the dissociative recombination. After Rydberg pathways were included in the theoretical de- 
scription, along with the Jahn-Teller coupling mechanism that excites the vibrational angular mo- 
mentum mode of the ion, DR theory was able to resolve the discrepancy. Theory and experiment 
for this fundamental chemical rearrangement process has now progressed to the point that some 
energy ranges can even be compared at the level of individual resonance features. Jungen and 
Pratt have independently demonstrated that the overall value of the DR rate coefficient can be ac- 
curately determined from a simplified model once the Jahn-Teller capture mechanism is included. 

Also in 2003, mid-infrared laser lines at wavelengths near 7 microns in laboratory hydro- 
gen/rare gas supersonic plasmas were observed at Berkeley.— Interestingly, strong IR emission 
from several massive star- forming regions is observed in a similar wavelength range of the spec- 
trum. Later, these laser lines in the Berkeley experiments were assigned to transitions between 
metastable H3 Rydberg states, as had been suggested by some detailed theoretical calculations.— 
A lasing mechanism was also proposed: the population inversion is generated by recombination of 
the ubiquitous molecular ion with low-energy electrons. Studies of flowing afterglow plasmas 
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by Glosik et al. suggest a three-body "collision assisted recombination" mechanism, rather than 
a simple two-body process because of the high (10^^ cm^) He gas density that is present in the 
supersonic discharge source.— 

More recently, experiments that study lasing in other energy ranges and in systems of other iso- 
topologues such as D3 in similar experimental conditions have been renewed. This has motivated 
us to extend our previous studies to this wavelength range at around 10-13 micron and to calculate 
the properties of lasing transitions between the Rydberg states of H3. An extension of our previous 
study to treat Rydberg states of the other isotopologue D3 is also presented. 



Our theoretical approach to the Rydberg states of H3 is based on multi-channel quantum defect 
theory (MQDT), one of the most successful techniques for treating Rydberg states in ab initio 
theory. This approach has been detailed in previous work,— so it will only be reviewed briefly 
here. 

In our studies, the model of studying molecular Rydberg energy levels of H3 treats the molecule 
as a Rydberg electron attached to the ion. The interactions between the Rydberg electron and 
the ion are described by body-frame quantum defects (or the equivalent reaction matrix elements^) 
that depend on the nuclear geometry. In the MQDT approach, a rovibrational transformation can 
be applied to construct the lab-frame ^-matrix using the body-frame quantum defect and the rovi- 
brational wave functions. For p-wave Rydberg states, the body-frame quantum defect parameters 
can be extracted from ab initio electronic potential surfaces. For higher orbital angular momentum 
states (Z > 1), a long-range multipole potential model is adopted. The rovibrational transformation 
can be formulated as follows: 



Here is an element of the laboratory-frame ^-matrix, which can be used to solve for eigenener- 



Method 




(1) 
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gies E of H3 by solving the following equation, which is the condition to kill exponentially growing 
components of the wavefunction at 0°: 

det|tan(;rv)+is:| =0. (2) 

The laboratory-frame eigenchannels \i) and the body-frame eigenchannels |a) are connected by 
the unitary transformation matrix Uia = {i\ oc), using the rovibrational wave functions of the 
ion core. To calculate these rovibrational wave functions, an accurate potential energy surface 
of is used,— 1^ and the three-body Schrodinger equation is solved within the hyperspherical 
adiabatic representation. In a recent paper, rovibrational energy levels of are calculated and 
compared with experiment with an accuracy at about 0.2 cm ^ Observe that Polyansky and Ten- 
nyson achieved an accuracy of 0.02 cm~^ using Jacobi coordinates. Their higher accuracy is due 
to the inclusion of nonadiabatic effects by using different effective reduced masses for vibration 
and rotation degree of freedom. Because the implementation of their procedure in hyperspherical 
coordinates is unclear, we have not attempted to reach this higher level of accuracy in the present 
calculations. However, the permutation symmetry of the rovibrational wave functions can be easily 
set up in hyperspherical coordinates, which is an important aspect of the rovibrational transforma- 
tion. Also, the accuracy of the computed Rydberg state energies of H3 is mainly limited by the 
accuracy of the body frame quantum defects, which yields uncertainties of typically a few cm^^ 
Therefore, the accuracy of the hyperspherical representation is adequate for our present purposes. 

The hyperspherical coordinates {R,0,(p} used in our approach are of the Smith-Whitten type,— 
which can be defined by the three interparticle distance ri2, and r^i through the relations: 

= 3-^/4i?[l + sin0sin((p- V6)]^/^ 

= 3"l/4i^[l+sin0sin((p-5;^/6)]^/^ (3) 
rgi = 3"i/4i?[l + sin0sin((p + ;r/2)]^/^ 

Together with the Euler angles a, /3 and 7, the three-body system can be described in the body- 
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frame. Similar to the usual Bom-Oppenheimer approximation, the adiabatic approach treats the 
hyperadius R initially as an adiabatic variable, and diagonalizes the Hamiltonian in all other de- 
grees of freedom (such as the hyperangles, Euler angles and spin degrees of freedom) yielding a 
set of adiabatic potentials and channel functions. The adiabatic corrections and couplings are later 
included using the "slow variable discretization" method.-^^^ 

One of the advantages of adopting this choice for the hyperspherical coordinates is that the 
basis functions used to discretize the Hamiltonian with the proper permutation symmetry can be 
easily constructed as, 



jm2K+ 



uj (0) - (-1)^"+^" 



SI 



(4) 



where uj (0) are a set of fifth-order basis splines which is unaffected by permutations. Here, the 
rotational part (OJ, /3, 7) is given by. 



2A^+ + 1 
8;r2 



(5) 



where D^j+j^+ are the Wigner D functions of the Euler angles. The phase of the Wigner func- 
tion is chosen as by Varshalovich et al.^^ is the total angular momentum of the ion, is 
the projection of onto the laboratory frame's z-axis, and m+ is the projection onto the body 
frame's Z-axis. 4>p is symmetry-adapted combinations of nuclear-spin functions for three spin half 
fermions defined as in a previous paper.— F — {A,E} represent the the symmetry representations, 
where gi = for T =A and gj = ±1 (ortho) for T = E (para). The permutation symmetries for the 
basis functions chosen for each degree of freedom are shown in[T] 

Under the condition that m2 + gi = 3n for even K^, and m2+ gi = 3n + 3/2 for odd K^, it is 
easy to show that the basis function obeys the permutation symmetry required for three identical 
fermions: 

(6a) 



'jm2K+ ' 
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Table 1 : Permutation symmetry for basis functions of different degrees of freedom. 



Permutation 
Operation 






4.r 

gl 


Pn 
P23 
P31 

PnP3i 
P12P23 


^i4n/3^-im2(p 
^i2n/3^-im2(p 
^i2n^-im2(p 
^i2K/3^im2(p 
^i4n/3 ^im2<p 


^K+m+ 


gi47tgf/3^r ^ 
^2ltgi/3^T ^ 
^ilKgi^Y ^ 
^i2Kgi/3^Y 



and 



jm2K^ jm2K^ ' ^ ^ 



where 

£^=1-Pl2-P23- P3l + P12P3I + Pl2P23- (V) 



Rydberg transitions of H3 in the 10 - 13 micron range 

The method described in last section has been applied to calculate 3p and 3d Rydberg states of 
H3, showing good agreement with experiments. The Ad ^ Ap and 6d 5p Rydberg transitions 
were used in Ref. to assign mid-infrared laser lines at wavelengths near 7 microns in laboratory 
hydrogen/rare gas supersonic plasmas. Here, the Rydberg transitions near 10-13 microns are 
calculated and shown m\T\ These transitions are mainly Id — t- 6p, 6d — > 6p and 5d — t- 6p Rydberg 
transitions. 



3/7;r Rydberg states of D3 

Using the method developed to calculate the Rydberg state energy levels for H3, we have also 
calculated energy levels for Spj^ Rydberg states of D3. The first step is again calculating the 
rovibrational states of the ion. In this calculation, the ionic potential surface for is adopted 
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Figure 1: Calculated nd — )■ n' p transitions of H3 Rydberg states at energy ranging from 750 to 1200 
cm~^ The y axis shows the theoretical Einstein B-coefficients in units of 10^^ (m/Js^). 
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from calculations done by Cencek et. al.,— »^ while the same quantum defects as in the case of 
are utilized for the Rydberg state calculation. This should be a good approximation since the 
quantum defects were calculated under the usual Bom-Oppenheimer approximation, where the 
masses of nucleus are assumed to be infinite. Nevertheless, the rovibrational energy levels of 
are calculated using nuclei mass of deuterium. 

A major difference between the calculations of rovibrational states of and is the dif- 
ferent permutational symmetries for the two species: the deuterium nuclei are bosons while the 
hydrogen nuclei are fermions. The symmetry-adapted combinations of nuclear-spin functions 4>^^ 
are constructed in the same way as given in by Kokoouline et. al. Here F = {Ai , A2, £} represents 
the symmetry representations of spin permutation group, and ^/ = for Ai (ortho) and A2 (para) 
symmetry, while ^/ = ±1 for E (meta) symmetry, since E representation is two dimensional. 
The permutation symmetry of these spin functions are tabulated in|2l 

Table 2: Permutation symmetry of symmetry-adapted spin functions for D^. 

Permutation 4>^/ (ortho) 4>g/ (para) 4>^^ (meta) 
Operation {gi = 0) {gi = 0) (^/ = ±1) 

~pj^ g'47rg//3^Ai _(,i47tgi/3f^A2 ^i47tgi/3f^E 

Si Si Si 

/'23 g'2?rg//3^^i _gi27tgi/3^A2 ^lugi/'i^E^ 

Si Si Si 

e'^^Si^^^^ _g'2;r^/^A2^ ^ilngj^E^^ 

PnP-il e'^^'SlP^^I ^i2Kg,/3^A2 ^i2ng,/3^E 

PnP23 g'4^g//34)^2 g'4;rg//34>£ 



The total nuclear-molecular function (including other degree of freedom such as rotation and 
vibration) should obey the permutation symmetry of three boson system, for example. 



Pn^''^i,^" = ^''XJ", (8) 



and 



^4>^"-> = 4>^"->, (9) 

]ni2K^ jm2K^ ' ^ ^ 



where ^ = 1 +^12 +^'23 +^'31 +/'i2/'3i + A2/'23- 
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Therefore, the basis functions are constructed as, 



4)" 



SI 



forr = Ai or E, and, 



(10) 



4)" 



N+m+Fg, 



(11) 



for r = A2, where g/, mi and satisfies m2+gi = 3n for even K^, and m2 + ^/ = 3n + 3/2 for 
odd K+. 

Using these numerical basis states having the appropriate permutation symmetry, the rovi- 
brational states of are calculated and compared with experimental results in|3l The r.m.s. 
difference between our calculation and experimental results is about 0. 1 1 cm ^ 

Using these accurate rovibrational states, a rovibrational frame transformation is applied to 
calculate the 3p;r Rydberg states of D3, and compared with experiment results— in HI From this 
table, the r.m.s. differences between experiment and our calculations are about 6 cm^ for almost 
all the results here. This might due to the quantum defect surface are optimal for H3, and the 
accuracy of our result might be improved by simply shifting the quantum defect by a small constant 
amount. 
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Table 3: Comparison of the calculated infrared transitions with the experimental values. 
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(cm-') 


{0,1,1} 


1 







{0,0,0} 
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1 
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0.086 
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1 


1 


£ 
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1 
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32.323 


1846.166 


1846.256 


0.090 
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1 


£ 
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1 





1979.123 


43.605 


1935.518 


1935.609 


0.091 
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2 
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£ 


{0,0,0} 


1 


1 


1967.997 
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1935.675 


1935.765 
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{0,1,1} 
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{0,0,0} 


2 





1887.976 


130.578 


1757.397 


1757.479 


0.082 
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1 


1 


£ 
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2 


1 


1878.488 


19.364 
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2 
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{0,1,1} 


2 





Ai 
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2 
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3 
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3 


A2 
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3 


3 
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159.859 


1720.403 
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3 


2 


£ 
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2 
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2047.093 


85.625 
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3 
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£ 
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2 
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2082.031 


119.364 
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3 
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£ 


{0,0,0} 


2 
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2068.516 
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{0,1,1} 


3 


1 


£ 


{0,0,0} 


2 


1 
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119.364 
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0.095 


{0,1,1} 
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Ai 


{0,0,0} 


2 





2115.066 


130.578 


1984.488 


1984.583 


0.095 


{0,1,1} 


3 


2 


£ 


{0,0,0} 


3 


2 


2068.516 


215.904 


1852.612 


1852.700 


0.088 


{0,1,1} 


3 


1 


£ 


{0,0,0} 


3 


1 


2103.337 


249.337 


1854.000 


1854.086 


0.086 


{0,1,1} 


3 


2 


£ 


{0,0,0} 


3 


2 


2047.093 


215.904 


1831.188 


1831.274 
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{0,1,1} 


3 


1 


£ 
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3 


1 
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3 


3 


A2 
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3 
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3 


3 
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3 


3 
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2 
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3 
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1 


2 
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1 


1 
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2 


4 


£ 
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1 


1 


3662.204 
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2 


3 


A2 
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2 


2 


£ 
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1 


1 
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1 


3 


Ai 
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1 


2 
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2 


1 
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1 


1 
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2 


2 


3718.205 
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4 
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2 


1 


3662.204 
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3542.840 
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{0,2,2} 


2 


3 


Ai 
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130.578 
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{0,2,2} 


2 


2 


£ 
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2 


1 
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3663.714 
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2 


1 


£ 
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2 


2 
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3720.984 
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2 


4 


£ 
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3 


1 
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